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ABSTRACT 



A two-component phenomenological model developed originally for £ Puppis is revised in order to model the outflows of 
late-type O dwarfs that exhibit the weak-wind phenomenon. With the theory's standard parameters for a generic weak- 
wind star, the ambient gas is heated to coronal temperatures w 3 x 10 6 K at radii J> 1AR, with cool radiatively-driven 
gas being then confined to dense clumps with filling factor w 0.02. Radiative driving ceases at radius w 2.1i? when the 
clumps are finally destroyed by heat conduction from the coronal gas. Thereafter, the outflow is a pure coronal wind, 
which cools and decelerates reaching oo with terminal velocity w 980 km s~ . 
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1. Introduction 

The X-ray emission from O stars (Harnden et al. 1979) is 
now generally agreed to arise from numerous shock fronts 
distributed throughout their winds. An early theory of such 
X-ray emitting winds (Lucy & White 1980; LW) was based 
on a two-component phenomenological model for the fi- 
nite amplitude state reached by unstable line-driven winds. 
Subsequently, the fundamental approach of computing the 
growth of the instability using the equations of radiation 
gas dynamics was pioneered by Owocki et al. (1988) and 
Feldmeier (1995), albeit with the then necessary restric- 
tions to 1-D flow and simplified radiative transfer. 

A question meriting further research is how and where 
this wind-shock model fails. According to LW, failure oc- 
curs at the low mass-loss rate ($) of a main sequence BO 
star, because the assumption of rapid radiative cooling of 
shocked ambient gas then breaks down for blob velocities 
<; 10 3 km s _1 , resulting in the heating of the blobs and 
consequent loss of line-driving. They conjecture that 'there- 
after, the relative motions of the two components dissipate 
and the smoothed wind coasts out to infinity' In effect, LW 
suggest that a wind that is initially radiatively driven con- 
verts into one that relies on thermal pressure to reach oo - 

1. e., a coronal wind. 

In addition to this question's intrinsic interest, it is no- 
table that the locus of this expected failure coincides with 
that of stars exhibiting the weak-wind phenomenon, (e.g., 
Marcolino et al. 2009; M09). Accordingly, this paper elab- 
orates LW's conjectures for the outflow from a weak-wind 
star. 

2. A Multi-zone wind: Zone 1 

The two-component model must be generalized to remove 
the assumption of instantaneous cooling of shocked gas and 
to incorporate blob destruction at finite radius. To achieve 
these aims, a multizone model is adopted, with each zone 
corresponding to different physical circumstances. 



Zone 1 starts just beyond the sonic point and ends when 
the isothermal-shock assumption is no longer justified. We 
assume instability has grown to full amplitide and adopt 
the LW description in which radiative-driven blobs (b) in- 
teract dynamically with a low density ambient (a) medium. 
Apart from the infinitesimally-thin cooling zones at shock 
fronts, both components are in thermal equilibrium with 
the photospheric radiation field, so that T a .b = T eq . 

2.1. Blobs 

The blobs can be identified with the clumps that are now 
a standard and spectroscopically-required feature of diag- 
nostic codes for O-star winds (e.g., Bouret et al. 2005). In 
such codes, the clumps are assumed to obey the /3— velocity 
law 
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where R is the photospheric radius and v^, the terminal 
velocity, is determined from the violet edges of P Cygni 
absorption troughs. 

Given the wide use of the /3-law, this now replaces LW's 
Eq. (10). But here, since Eq. (1) ceases to apply when r > 
r$, the blobs' destruction radius, is not an observable. 
The highest velocity at which UV absorption is detected is 
a measure not of Voo but of Vb(rs). 

For given /3, diagnostic modellers choose the clumps' 
filling factor fa and mass-loss rate $b so that absorp- 
tion troughs have their observed strengths. For a strong 
line at frequency i>q, this typically requires that, despite 
dumpiness, a continuum photon emitted between vq and 
i^o (1 + Vqo/c) has small probability of escaping to oo, and 
so most of the photon momentum in this interval is trans- 
ferred to the clumps. In an LW wind, essentially the same 
requirement arises as a consistency criterion: since the am- 
bient gas is assumed not to be radiatively driven, it must be 
shadowed by the blobs. The optical depth criterion adopted 
by LW is that n(r) > 1.5 for all r, where T\ is given in 
Eq.(ll) of LW. For the (3— velocity law, the function x m 
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their n formula becomes 

X (x; P) = l(l~ x) 1 - 2 ' 3 /W + x-{l + f3)x 2 } (2) 
where x = R/r. 



2.2. Dynamics 

In zone 1, the dynamical interaction of the blob and am- 
bient components is treated exactly as in LW. For more 
recent treatments and applications, see Howk ct al. (2000) 
and Guo (2010). 

The equation of motion obeyed by the blobs is 



v b~r = 9R~ 9d - 9 
dr 



(3) 



where gn,gD, and g are the forces per gram due to radi- 
ation, drag, and gravity, respectively. The drag force go 
retards the blobs but accelerates the ambient gas and is 
the means by which photon momentum is transferred to 
this component. The resulting equation of motion of the 
ambient gas is 



v — 

a dr 



1 dP a p b 
• — -j— + — 9D- 9 

Pa dr p a 



(4) 



Here p a ,b are the smoothed densities, and P a = a 2 p a with 
al = kT a /pm H . 

With the LW assumption of no mass exchange between 
the components, the two equations of continuity integrate 
to give 

$ 0i 6 = 4?rr 2 p a ,bV a ,b (5) 

where & a ,b are constants whose sum is the star's mass-loss 
rate <I>. 

The drag force mgD on a blob of mass m is computed 
using De Young and Axford's (1967) theory of inertially- 
confined plasma clouds - see Sect. II b) in LW. The resulting 
formula is 



C D p a U 2 A b 



(6) 



where Ab — nv 2 is the blob's mean cross section, U = 
Vb — v a is the blob's velocity relative to the ambient gas, 
and the drag coefficient Co — 1-519. 

2.3. Filling factors 

At a point (vb,v a ,r) in an outward integration with spec- 
ified $ a ,f>! the smoothed densities p a b are given by Eq. 
(5). The ambient density is then p a — p a /fa, where f a 
is the filling factor of the ambient gas. Correspondingly, 
the mean density of the stratified De Young - Axford blob 
is pb = pb/ h- Now, in the absence of a void component, 
fa + fb = 1, and so only one of f a and /& is indepen- 
dent. To determine fb, say, we must iterate. Solution by 
repeated bisection is adopted, starting with upper and 
lower limits fu = 1 and fu = 0. Then, with the esti- 
mate fb = (fu + /l)/2, the blob's volume Vb is computed 
from LWs Eq.(5). The mean density of the blob is then 
p b = ra/Vb, corresponding to f b = p b /p b . If fb < fb, the 
new upper limit is fu — f b - On the other hand, if fb> fb, 
the new lower limit is fu — fb- The iterations continue until 
fu ~ fu < 10~ 7 . Then, with the resulting converged value 
of fb, all quantities required to continue the integration can 
be evaluated. 



2.4. Switch criterion 

The assumption of instantaneous cooling breaks down at 
low densities because the cooling rate per unit volume 
C oc p 2 . If the cooling time scale t c increases to the ex- 
tent that a parcel of shock-heated gas encounters another 
shock before cooling back to T eq , then shock-heating raises 
the mean temperature of the ambient medium. An approx- 
imate criterion for this transition to zone 2 is derived as 
follows: 

First, consider radiatively-cooled flow of monatomic gas 
(7 = 5/3) emerging from a steady shock (see Fig.l in Draine 
& McKee 1993). Since this flow is subsonic, the pressure 
gradient may be neglected in comparison to that of tem- 
perature. Thus, in the shock's frame, 



dlnT 
dr 
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(7) 



The cooling timescale is therefore t c — £ c /v — 5/2 x P/C. 

Now consider flow into the bow shocks. The entire mass 
p a of ambient gas in unit volume is shocked in time interval 
U = Pa/Nbjb, where Mb = fb/Vb is the number density of 
blobs, and jb ~ p a UAb is the mass flow rate through each 
bow shock. Hence U w 4/3./), x a/U. 

The criterion for switching from zone 1 to zone 2 is then 
simply t c >t l . 



3. A Multi-zone wind: Zone 2 

The outward integration of the wind continues in zone 2 
with the same basic model except that T a b 7^ T eq . The 
ambient gas is now heated by being repeatedly shocked, 
and the blobs in turn gain heat by conduction from the 
ambient gas. Zone 2 ends when the blobs can no longer 
achieve thermal equilibrium. 



3.1. Blob survival 

The survival of blobs (clumps) in stellar winds has similar- 
ities to that of clouds in the interstellar medium. In that 
context, Cowie & McKee (1977) studied the evaporation 
of a spherical cloud embedded in a hot tenuous medium. 
Importantly, they treated the saturation of heat conduc- 
tion when the electron mean free path in the surround- 
ing medium is <; the cloud's radius and estimated, under 
the assumption of steady outflow, the reduced evaporation 
rate. In a companion paper (McKee & Cowie 1977; see also 
Graham & Langer 1973), they consider the effects of ra- 
diative losses, finding that evaporation is replaced by con- 
densation if the losses exceed the heat input from the hot 
gas. However, numerical calculations by Vieser & Hensler 
(2008) cast doubt on the assumption of steady outflow. In 
the case of saturated conduction, they find a further reduc- 
tion of evaporation rate by a factor ~ 40 due to changes of 
the cloud's environment caused by the outflow. 

Given the evident difficulty of reliably predicting when 
cool gas is eliminated by its interaction with surrounding 
hot gas, a simple prescriptive approach is adopted here: As 
in zone 1, the blobs retain their fixed mass m throughout 
zone 2. However, when the heat input from the ambient gas 
exceeds their maximum cooling rate, the blobs are assumed 
to merge instantly with the ambient gas. 
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3.2. Heating and cooling of blobs 

In zone 2, the blobs are surrounded by shock- heated gas and 
so will be heated by thermal conduction. But if the ambi- 
ent gas reaches coronal temperatures, heat conduction is 
flux-limited. Moreover, conductivity may be suppressed by 
magnetic fields. An approximate formula interpolating be- 
tween the classical and saturated limits and incorporating 
a suppression factor (f> is derived in Appendix A. 

If Li n is the rate of heat flow from the ambient gas, a 
blob will achieve thermal equilibrium at T b > T eq if the 
enhanced radiative cooling rate 

AC b = {n e n H ) b (A(T b ) - A(T eq )) x V b = C m (8) 

Here A(T) is the optically-thin cooling function, and the 
blob is treated as isothermal and of uniform density. (But 
note that LW's definition of p b is such that A£ b is exact for 
the density stratification of an isothermal De Young- Axford 
blob.) 

Because A(T) reaches a maximum at Tf(K) = 5.35 dex 
(Dere et al. 2009), the solution of Eq.(8) with T b < T t is 
appropriate as the blobs are heated to above T eq . When 
T b reaches 7^, the corresponding £ in is the maximum 
value consistent with thermal equilibrium. Any further in- 
crease in Ci n cannot be matched by increased cooling. 
Accordingly, we take this as the point beyond which the 
blobs cannot survive. 

Note that if conduction is completely suppressed (<fr = 
0), then d n = and the solution of Eq.(8) is T b = T eq . 
The blobs therefore survive, and zone 2 extends to oo. 

3.3. Heating and cooling of ambient gas 

According to LW, the rate at which energy is being dissi- 
pated per unit volume is 

Q ~ 9d p b U (9) 
Dividing by J\f b , we find that the rate per blob is 

Qb = \p a U 2 xUx C D A b (10) 

showing that in unit time each blob's bow shock dissipates 
the kinetic energy in a column of inflowing gas of length 
U and cross section CdAi>. In zone 1, this dissipated en- 
ergy is radiated by a thin cooling layer, and so Q deter- 
mines the sum of these layers' frequency-integrated emis- 
sivities - Eq.(7) in LW. But in zone 2 where t c > t i} we 
jump to the opposite limit, treating dissipation as a heat 
source distributed uniformly throughout the ambient gas, 
and similarly for cooling. Accordingly, the energy equation 
for stationary flow of the monatomic ambient gas is 

(dP a 5 2 dp a \ 2 • • 

Note that since f a w 1 terms arising from radial changes in 
f a have been neglected. 

The total cooling rate per unit volume of ambient gas, 
Ctot, is the sum of the losses due to radiative cooling and 
to conduction into the blobs. Thus, 

Ctot = (n e n H ) a A(T a ) + Af b £ in (12) 



Integration of Eqs. (4), (5) and (11) continues until T b = Tf, 
at which point the blobs are deemed to merge instantly with 
the ambient gas (Sect. 3.1). Accordingly, the transition from 
zone 2 to zone 3 occurs at S, a surface of discontinuity (e.g., 
Landau & Lifshitz 1959), across which the fluxes of mass 

J = J a + Jb = faPaVa + IbPb^b (13) 

momentum 

n = f a (Pa + PoVl) + fb(Pb + Pbvl) (14) 

and energy 

F = Ja Qt£ + \a 2 }j + J b {^v 2 b + |ag) (15) 

are continuous. Zone 2 thus ends at r$ with the evaluation 
of J, n and F. 

4. A Multi-zone wind: Zone 3 

The outward integration continues in zone 3, but now the 
blobs have disappeared, leaving a single fluid component 
with no driving force (go = 0) and no heat input (Q = 0). 
The initial conditions for the resulting ODE's are obtained 
from the continuity across S of J, LI and F. Thus, v, p, and 
T for the flow emerging from S + are given by 

pv = J (16) 

P + pv 2 = LI (17) 

and 

jQf + ty-F (18) 
where J, LI, and F are given by Eqs. (13)-(15). 

4.1. Solution branches 

From Eqs. (16)-(18), we readily derive the quadratic equa- 
tion 

v 2 - 2u p v + u 2 e = (19) 

where u p = 5LI/8J and u e = y/(F/2J). The two solutions 
are 

v± = Up ± \Jv%~-u% (20) 

The corresponding temperatures T± are derived from the 
isothermal sound speeds given by 

2 kT± 2 2 

a± = = u e - -u p v± (21) 

pm H 5 

and the densities are p± = J/v±. 

If Up — u e , the two solutions coincide. When this hap- 
pens, v = u p = y/(5/3) a - i.e., the outflow at S + is exactly 
sonic. If Up > u e , the solutions are real and distinct. The 
v + solution is supersonic (M + branch) , and the V- solution 
is subsonic (M_ branch). Note that the M+ branch has a 
singularity at v p /v e = 5/4, at which point a + = 0. Mach 
numbers for the two branches are plotted against v p /v e in 
Fig. 1. 

The M branch corresponds to S being the locus not 
only of merging but also of a stationary shock front. This 
branch would perhaps be appropriate if there were a pre- 
existing slower wind (cf. Macfarlane & Cassinclli 1989), but 
this is not the circumstance evisaged here. Instead, there- 
fore, the M + branch is selected since this corresponds to a 
high speed two-component flow at .S* - emerging at S + as a 
single-component supersonic flow. 
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Fig. 1. Mach number M as a function of u p /u e along the su- 
personic (-M+) and subsonic (M_) solution branches. When 
w p /ti e = 5/4, M+ = co and M_ = 1/^5 = 0.447. 



4.2. Dissipation at S 

In addition to the roots of Eq.(19) being real, a further con- 
dition is mandatory: the transition from S~ to S + must be 
such that kinetic energy is dissipated (entropy production) 
and not the reverse. For the two branches, kinetic energy 
is thcrmalized at the rates 



1 



L S = o(*a<+«-$«±) 



(22) 



whose positivity must be checked. 

Note that the kinetic energy dissipated at S is not ra- 
diated away by a thin cooling zone. Instead, this energy 
contributes to the flow's enthalpy at S + , which then does 
PdV work in the subsequent expansion. 

4.3. Outward integration 

The solution for the single-component gas in zone 3 is ob- 
tained by integrating the equations of motion 



dv 



IdP 



"T = -I 9 

dr p dr 



continuity 
and energy 



4:irr 2 pv 



dP 

dr 



5 2 dp 
3° * 



n e n H A(T) 



(23) 
(24) 

(25) 



The initial conditions at r$ are v + , p + and T + derived in 
Sect.4.1. 

This integration continues to r = oo. However, this is 
only possible if the energy density at S + is sufficient to over- 
come both the remaining potential barrier and the cooling 
losses. If not, a stationary, spherically-symmetric wind so- 
lution of this type does not exist. 



5. An example 

To illustrate the ideas presented in Sects. 2-4, the solution 
for a generic weak-wind star is now described in detail. 



5.1. Standard parameters 

The model has several parameters, for which standard val- 
ues are now adopted. Given their uncertainty, sensitivity to 
changes are reported in Sect. 6. 

Because the theory does not predict this is derived 
from previously-tabulated mass fluxes (Lucy 2010b; LlOb). 
The chosen model has T c « = 32.5kK and logg = 3.75, 
consistent with the weak- wind stars £ Oph and HD 216532 
- see Table 3 in M09. The model's mass flux J(gm cm~ 2 
s- 1 ) = -7.11 dex. 

The star's mass M = 2A.1M.Q is determined by find- 
ing the point on the ZAMS from which the evolution- 
ary track during core H-burning has log g = 3.75 when 
T cS = 32.5kK. This point is reached after 5.75 x 10 6 yrs 
when R = 10.83i? Q and the luminosity L = 1.18x 10 5 Lq = 
4.52 x 10 38 erg s _1 . The assumed composition is X = 
0.70, Z = 0.02. 

With R and L determined, $ = 4ttR 2 J = 8.80 x 
10~ 9 A4q yr -1 = l.lOL/c 2 . This theoretical <& derives from 
the constraint of regularity at the sonic point (v = a) in 
the theory of moving reversing layers. In the weak-wind 
domain, this theory's predictions exceed the highly uncer- 
tain (±0.7 dex) observational estimates of M09 by w 0.8 
dex but are lower than the Vink et al. (2000) formula by 
w 1.4 dex (Lucy 2010a; LlOa). 

For the parameters in Eq.(l), we adopt the 
observationally-supported O-star values f3 = 1 and 
^oo = 2.6iw(i?) = 2394 km s _1 . 

The mass m of the blobs must also be specified. Recent 
modelling of O-star spectra finds that 'In most cases, 
clumping must start deep in the wind, just above the sonic 
point' (Bouret et al 2008). We therefore retain LW's as- 
sumption that blobs form at or near the sonic point and 
have diameters comparable to H pi the local scale height. 
At v = a in model i325#375, p = 4.27 x 
and H n 



10 14 gm cm 3 



9.86 x 10 8 cm, so that the crude LW estimate is 
gm. 



m = 2 x 10 13 _ 

The ratio ry = < ! ) &/ < E > must also be specified. Following 
LW, we determine r\ by imposing the constraint that r m = 
1.5, where 

r m — min [n (r)] in zone 1 (26) 

Typically, the minimum occurs at the end of zone 1 
where inertial confinement is greatest. In zone 2, shadowing 
rapidly becomes irrelevant since the rapid rise of T a - see 
Fig. 3 - destroys driving ions. 

Finally, the conductivity suppression factor <fi intro- 
duced in Appendix A must be specified. As standard value, 
we set <j) = —1.0 dex, a moderate degree of suppression 
compared to estimates for galaxy clusters (e.g., Ettori & 
Fabian 2000). 

5.2. Zone 1 

In this high-density zone close to the photosphere, both 
components are assumed to be in thermal equilibrium with 
the star's radiation field, a condition approximated by set- 
ting T eq = 0.75T e fj, as in L10a,b. 
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With the assumptions of isothermal flow, specified Vb, 
and no mass exchange between components, the structure 
of zone 1 is obtained by integrating the ODE 



dllLV a 

d\nr 



2< , Pb 
1- — 9d 

r Pa 



(27) 



The outward integration starts, as in LW, with Vb = 150 
km s _1 and v a — 100 km s _1 , a point sufficiently beyond 
the presumed onset of dumpiness that the two-component 
state may be regarded as established. The starting radius 
from Eq.'(l) is n = 1.067i?. 

Eq.(27) has a singularity when v a = a a . Since the in- 
tegration starts with v a > a a , this singularity only arises 
if insufficient drag gr> causes the flow to decelerate. A pa- 
rameter set for which this happens does not admit a steady 
wind of this type. 

As shown in Fig. 2, the standard parameters result in an 
outflow of ambient gas that accelerates throughout zone 1 . 
This continues until the switch to zone 2 is triggered by the 
onset of the inequality t c > t% - see Sect. 2. 4. This occurs at 
r/R = 1.28, with v b '= 528 km s" 1 and v a = 325 km s" 1 . 
The relevant time-scales are t c = t{ = 2.0 x 10 3 s, which 
are <§C the local flow time-scale, r/vb = 1.8 x 10 4 s. 

The post-shock cooling rate C required in calculating t c 
is given by n e n#A(T), where A(T) is the optically- thin 
cooling function for photospheric abundances tabulated by 
Dere et al. (2009). This rate is computed at the apex of the 
bow shock with n e = 1.18njj, corresponding to complete 
electron-stripping. 

At the end of zone 1, the post-shock temperature has 
risen to 6.0 x 10 5 K, so X-ray emission from zone 1 is neg- 
ligible. 

5.3. Zone 2 

With the isothermal assumption dropped, the structure of 
zone 2 is determined by Eqs. (4) and (11). With dependent 
variables v a and T a , the ODE's to be integrated are 



{vl - al) 



din v a 
dlnr 



2 d In T a 



"2a 2 a + r ( j-g D - 9 



(28) 



and 



2 d\nv a d\nT a 

3 d In r d In r 



4 2 r 
3 + 3P~^ 



(Q-C tot ) (29) 



Since all variables are continuous at this transition, the in- 
tegration starts at the point (v a ,Vb,T a ,Tb 7 r) reached by 
the zone-1 integration. 

Eqs. (28) and (29) are a pair of algebraic equations for 
the two derivatives. The determinant of the coefficients' 
matrix is zero when v a = v / (5/3)a a - i.e., at the adiabatic 
sonic point. If this singularity is ecountered, the parameters 
are inconsistent with the conjectured wind structure. 

Fig. 2 shows that, with the standard parameters, the flow 
continues to accelerate throughout zone 2 reaching v a = 940 
km s _1 at rs = 2.14_R, at which point Vb = 1277 km s _1 . 

The corresponding temperature structure predicted for 
zone 2 is shown in Fig.3. At the start, T a f, = T eq = 24.4kK. 
Thereafter, shock-heating of the ambient component over- 
comes radiative, conductive and adiabatic cooling to give a 
rapidly increasing T a , reaching the coronal value 10 6 K at 
r = IMR and 3.7 x 10 6 K at r s . 
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Fig. 2. Velocities of blobs (vb) and ambient gas (v a ) as func- 
tions of radius. Zone boundaries are indicated. The surface of 
discontinuity S where blobs merge with ambient gas occurs at 
r/R = 2.14. 



The profile for T b shows discontinuous jumps at the be- 
ginning and end of zone 2. These result from non-monotonic 
variations of A(T). For example, A's peak at T^(K) = 5.35 
dex is preceded by lower peak at 5.00 dex. Accordingly, 
after reaching Tb(K) = 5.00 dex, a slight increase in £ in 
results in a discontinuous jump to Tb(K) = 5.18 dex, fol- 
lowed quickly by blob destruction when Tb = Tf. Because 
of these jumps, the radiative driving of the blobs, which 
is ultimately reponsible for T a 's increase to coronal values, 
occurs mostly between T& = 40 and 90kK. 

Blob temperatures are derived algebraically from Eq. (8) 
on the assumption that blobs adjust instantaneously 
to thermal equilibrium. At rs, the heating time scale 
1.5 nkTf x Vb/Ci n = 0.9 x 10 2 s compared to the flow time 
scale r/vb = 1.3 x 10 4 s. 

In computing cooling rates for blobs, we set n e = 
1.12n#, corresponding to metals being stripped of ~ 2 — 3 
electrons. 



5.4. Surface of discontinuity S 

At S~, the blobs have filling factor /& = 0.024, velocity 
v b = 1282km s" 1 and temperature T b = 2.24 x 10 5 K. 
The corresponding values for the ambient component are 
f a = 0.976, v a = 944km s^ 1 and T a = 3.66 x 10 6 K. After 
merging, the flow at S + has two possible solutions (Sect. 
4.1). For the rejected M_ solution, the flow emerges with 
v_ = 313km s~\ T_ = 1.90 x 10 7 K, corresponding to Mach 
0.48, and the implied rate at which kinetic energy is dissi- 

- 1 or 6.9 x 10" 6 L. 



pated L s = 0.31 x 10 34 erg s" 

For the selected M + solution, the flow emerges with 
v + = 1094km s-\ T+ = 2.58 x 10 6 K, corresponding to 
Mach 4.6, and the implied dissipation rate Lg = 0.92 x 10 32 
erg s _1 or 2.0 x 10~ 7 L. 

Notice that v + € (v a ,v b ), as expected if 5* is the locus 
only of merging. In contrast, u_ < v a , so there is a coinci- 



5 



Lucy: Coronal winds 



l 


t 








6.5 






a . 


6 




\ 


zone 3 


5.5 


zone 1 


\ zons 2 




5 








4.5 




Temperature profiles 











0.1 0.2 0.3 0.4 0.5 

Log r/R 



Fig. 3. Temperatures of blobs (6) and ambient gas (a) as func- 
tions of radius. Zone boundaries are indicated. 

dent shock, as also indicated by the far greater dissipation 
rate Lg. 

5.5. Zone 3 

The single component flow emerging from S is a pure coro- 
nal wind: the only outward force is the gradient of thermal 
pressure. 

The structure of zone 3 is obtained by continuing the 
integration of Eqs. (28) and (29), but now with go = 0, Q = 
and Ctot — n e njjA(T). The initial conditions at r$ are 
v + and T + given in Sect. 5.4. 

A short segment of this outflow is plotted in Figs. 2 
and 3, showing that the flow decelerates and (inevitably) 
cools. For these standard parameters, the energy density at 
S suffices to overcome cooling and power escape to oo. At 
rf/R = 100, the flow has slowed to 984km s _1 , way beyond 
the local v esc — 92km s _1 

The temperature drops below the coronal value 10 6 K at 
r/R = 4.24 and to 10 5 K at r/R = 13.2. 

5.6. Emission measure 

With standard parameters, our generic weak-wind star is 
predicted to have a corona (T > 10 6 ) that extends from 
r\ = 1.35i? to r2 = 4.24i? and so will be an X-ray emitter. 
As a crude guide to detectability, we compute the emission 
measure of coronal gas 

/T2 

e = 4-7T / n e riH r 2 dr (30) 

Jri 

and its hardness parameter 

< kT > = Airs' 1 / kTn e n H r 2 dr (31) 
J n 

The results are s(c m - 3 ) = 53.51 dex and < kT >= 0.20 
keV. 



5. 7. Energy budget 

The global energy budget of this multi-zone wind is of in- 
terest. The input is the rate of working in zones 1 and 2 of 
gn, the force per unit mass acting on the blobs. This rate 
L W rk = 5.4 x 10 33 erg s _1 . 

The balancing output is Lm + Lw, where Lm is the rate 
at which matter gains kinetic and potential energy, and L\y 
is the wind's radiative luminosity. For the interval (r i7 r/), 
L M = 4.8 x 10 33 erg s _1 or 88.5% of L wrk . The remaining 
11.5% is accounted for by Lw, which comprises radiative 
losses from shock fronts in zone 1, cooling radiation from 
blobs and ambient gas in zone 2, and cooling radiation from 
the coronal flow in zone 3. 

For an idealized line-driven wind in which gas remains 
(by assumption) at T eq , PdV work is negligble so that 
Lm — L wr k- In contrast, for a pure coronal wind, L wrk = 0, 
so that Lm is entirely due to the PdV work of the hot gas. 
The relative contributions of these two mechanisms in this 
hybrid case is of interest. 

In answering this, we must first integrate Q from Eq.(9) 
over zones 1 and 2 to obtain the total dissipation rate L n = 
1.1 x 10 33 erg s _1 . The quantity L wrk -L D = 4.3 x 10 33 erg 
s _1 is then the contribution to Lm due directly to radiative 
driving. On the other hand, the contribution of PdV work 
by hot gas is L D — L w = 0.4 x 10 33 erg s _1 . 

A measure of the proximity of a hybrid- to a pure coro- 
nal wind is the ratio 

9 = (L D - L W )/L M (32) 

which = for a conventional line-driven wind and = 1 
for a coronal wind. With standard parameters, the multi- 
zone wind has 9 = 0.08, so direct radiative driving still 
dominates in accounting for Lm- 

A further quantity of interest is the integrated cooling 
rate of gas with T e > 10 6 K, since this is approximately 
the wind's X-ray luminosity. For zones 2 and 3, this gives 
L x w 3.4 x 10 31 erg s~\ so that L x /L w 0.76 x 10~ 7 L, 
similar to the ratio found for early-type O stars. 

6. Non-standard parameters 

The theory developed in Sects. 2-4 has several param- 
eters, each of which would either be predicted or ren- 
dered unnecessary if calculations could be carried out from 
first principles. Sensitivity of the results to these currently 
unavoidable parameters must therefore be investigated. 
Accordingly, sequences of solutions are now reported in 
which a single parameter is varied while keeping others at 
the standard values of Sect. 5. 

Key properties of the models are given in Table 1. The 
quantities reported are as follows: 

Col. 1: Sequence identifier. 

Col. 2: Exponent in Eq.(l), the velocity law. 

Col. 3: Log of total mass-loss rate in A4q yr _1 . 

Col. 4- Log of blobs' mass in gm. 

Col. 5: Log of conductivity suppression factor - see Eq. 
(A.9). 
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Fig. 4. Sensitivity of coronal winds to 0, the magnetic suppres- 
sion factor - see Eq.(A.9). Values of log <j> are shown. The vertical 
segments are the surfaces of discontinuity S. 



Col. 6: Fraction of mass- loss in blobs = $&/$ . 

Col. 7: Shadowing optical depth - see Eq.(26). 

Col. 8: Vb in km s _1 at the destruction radius r$. 

Col. 9: Maximum ambient gas temperate in 10 6 K. 

Col. 10: Log of emission measure in cm -3 - see Eq.(30). 

Col. 11: Hardness parameter in keV - see Eq.(31). 



6.1. Sequence I 

In this sequence, the conductivity suppression factor varies 
from 4> = 0.001, an extreme value but with observational 
support for galaxy clusters (Ettori & Fabian 2000), to <p = 
1, the value for a non-magnetized plasma. 

Not surprisingly, the predictions are highly sensitive to 
4>, and this might eventually be exploited diagnostically. 
With = 1, conductive heating of the blobs destroys them 
already at r$ = 1.57R where Vb = 873km s _1 . However, 
with 4> = 0.001, blobs survive out to r$ = 22. OR where 
v b = 2286km s" 1 . 

Table 1 also shows that coronal temperature and the 
hardness parameter increase as (f> — > 0. However, the emis- 
sion measure e remains <~ 53.4 — 53.5 dex after an initial 
sharp rise from 53.18 dex for 0=1. The predicted tem- 
perature profiles of the coronae as (j) varies are plotted in 
Fig.4. 

This sequence demonstrates the diagnostic potential of 
UV and X-ray data in constraining magnetic suppression 
of conductivity. The UV data measures the highest velocity 
at which wind matter transfers photon momentum to the 
gas and the X-ray data measures the hardness of coronal 
emission. 



A further diagnostic test provided by P Cygni lines is 
the weakness of emission components. As r$ decreases with 
increasing (j), the fraction of scattered photons occulted 
by the star increases and the emission component weak- 
ens. This effect was invoked for r Sco by LW in arguing 
that 'outflowing gas loses its ability to scatter UV radia- 
tion while still close to the star's surface.' Note that the 
weak-wind stars investigated in M09 all have C iv reso- 
nance doublets with weak or absent emission components. 
Diagnostic modelling of these stars would improve if UV 
scattering were truncated at finite radius. 

6.2. Sequence II 

As noted in Sect. 5.1, the $ of weak-wind stars is poorly 
determined. This sequence explores sensitivity to this un- 
certain parameter. 

When $ is increased above the standard value from 
LlOb, the blobs survive to higher velocities, and the higher 
coronal densities give the approximate scaling law s cx 
Interestingly, the quantities T max and < kT > are insensi- 
tive to $. 

The attempt to continue this sequence to lower $'s 
failed at —8.36 dex because the singularity in zone 2 dis- 
cussed in Sect. 5. 3 is encountered. This arises as follows: the 
sharp initial rise of T a in zone 2 causes M to decrease de- 
spite increasing v a - see Figs. 2 and 3. But as T a levels off 
M reaches a minimum and then rises again. Sequence II 
terminates for $ between —8.36 and —8.26 dex when this 
minimum falls to M = 1. For the solution plotted in Figs. 
2 and 3, this zone-2 minimum is M = 1.88 at r = 1A7R. 

6.3. Sequences lll-V 

In sequence III, the velocity-law exponent varies from ft = 
0.5 - rapid acceleration - to ft — 2.5 - slow acceleration. 
The standard value ft = 1.0 is approximately a stationary 
point as regards the coronal properties e and < kT >, so 
these are insensitive to ft. However, Vb(rs) is moderately 
sensitive. 

In sequence IV, solution sensitivity to the highly un- 
certain blob mass is explored. Fortunately, coronal prop- 
erties are only moderately sensitive, with e oc to -0 22 and 
< kT > cx to 03 . In regard to blob destruction, this occurs 
as expected at low velocities for small m. In consequence, 
sequence IV terminates for m{gm) between 11.8 and 11.9 
dex because the outflow in zone 3 is then unable to reach 
oo on account of negative energy density - see Sect. 4. 3. 

Finally, sensitivity to the shadowing parameter r m is 
investigated with sequence V. Again only moderate sensi- 
tivity is found. 

7. Conclusion 

The aim of this paper has been to investigate the struc- 
tural changes of O-star winds when $ decreases to the ex- 
tent found for the weak-wind stars. To this end, the two- 
component phenomenological model developed originally 
for ( Puppis is modified to incorporate LW's conjectures 
following the breakdown of that model's assumptions for 
r Sco. When applied to a generic weak-wind star, the re- 
vised model predicts that shock-heating of the ambient gas 
gives rise to coronal temperatures, that conductive-heating 
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Table 1. Solutions with non-standard parameters. 





a 
P 


log 


log vet 


log <f) 


V 


7"m 


Vb{'S ) 


T 

J- max 


log e 




I 


1.0 


-8.06 


13.3 


-3.0 


0.45 


1.5 


2286 


8.4 


53.43 


0.40 




1.0 


-8.06 


13.3 


-2.5 


0.45 


1.5 


2180 


8.2 


53.44 


0.39 




1.0 


-8.06 


13.3 


-2.0 


0.45 


1.5 


1991 


7.7 


53.44 


0.36 




1.0 


-8.06 


13.3 


-1.5 


0.45 


1.5 


1678 


6.2 


53.47 


0.30 




1.0 


-8.06 


13.3 


-1.0 


0.45 


1.5 


1277 


3.7 


53.51 


0.20 




1.0 


-8.06 


13.3 


-0.5 


0.45 


1.5 


1015 


2.1 


53.47 


0.13 




1.0 


-8.06 


13.3 


0.0 


0.45 


1.5 


873 


1.4 


53.18 


0.11 


II 


1.0 


-8.26 


13.3 


-1.0 


0.47 


1.5 


1109 


3.5 


53.27 


0.19 




1.0 


-8.06 


13.3 


-1.0 


0.45 


1.5 


1277 


3.7 


53.51 


0.20 




1.0 


-7.76 


13.3 


-1.0 


0.41 


1.5 


1509 


3.7 


53.92 


0.21 




1.0 


-7.46 


13.3 


-1.0 


0.37 


1.5 


1719 


3.6 


54.31 


0.21 




1.0 


-7.16 


13.3 


-1.0 


0.32 


1.5 


1918 


3.3 


54.70 


0.21 




1.0 


-6.86 


13.3 


-1.0 


0.27 


1.5 


2121 


2.9 


55.06 


0.19 


III 


0.5 


-8.06 


13.3 


-1.0 


0.68 


1.5 


1348 


3.6 


53.14 


0.17 




0.6 


-8.06 


13.3 


-1.0 


0.59 


1.5 


1380 


4.0 


53.37 


0.20 




0.8 


-8.06 


13.3 


-1.0 


0.50 


1.5 


1341 


3.9 


53.49 


0.20 




1.0 


-8.06 


13.3 


-1.0 


0.45 


1.5 


1277 


3.7 


53.51 


0.20 




1.5 


-8.06 


13.3 


-1.0 


0.37 


1.5 


1135 


3.2 


53.50 


0.17 




2.0 


-8.06 


13.3 


-1.0 


0.32 


1.5 


1033 


2.8 


53.45 


0.15 




2.5 


-8.06 


13.3 


-1.0 


0.29 


1.5 


962 


2.5 


53.40 


0.14 


IV 


1.0 


-8.06 


11.9 


-1.0 


0.20 


1.5 


889 


2.6 


53.81 


0.14 




1.0 


-8.06 


12.3 


-1.0 


0.26 


1.5 


989 


3.0 


53.74 


0.16 




1.0 


-8.06 


12.8 


-1.0 


0.35 


1.5 


1135 


3.4 


53.63 


0.18 




1.0 


-8.06 


13.3 


-1.0 


0.45 


1.5 


1277 


3.7 


53.51 


0.20 




1.0 


-8.06 


13.8 


-1.0 


0.55 


1.5 


1400 


4.0 


53.38 


0.21 




1.0 


-8.06 


14.3 


-1.0 


0.66 


1.5 


1503 


4.2 


53.23 


0.22 




1.0 


-8.06 


14.6 


-1.0 


0.72 


1.5 


1555 


4.4 


53.12 


0.22 


V 


1.0 


-8.06 


13.3 


-1.0 


0.41 


1.0 


1322 


4.2 


53.67 


0.24 




1.0 


-8.06 


13.3 


-1.0 


0.45 


1.5 


1277 


3.7 


53.51 


0.20 




1.0 


-8.06 


13.3 


-1.0 


0.48 


2.0 


1246 


3.3 


53.40 


0.17 




1.0 


-8.06 


13.3 


-1.0 


0.50 


2.5 


1225 


2.9 


53.31 


0.16 



eventually destroys the blobs, and that the resulting single- 
component flow coasts to oo as a pure coronal wind. Thus, 
in broad outline, the volumetric roles of hot and cool gas 
in O-star winds are reversed. In the now standard picture 
for a star such as £ Puppis the X-ray emitting gas occupies 
a tiny fraction of the wind's volume, with the bulk of the 
volume being highly-clumped cool gas with T ~ T eq . In con- 
trast, in the picture suggested here for the weak-wind stars, 
X-ray emitting gas fills most of the volume for r <; 1.3R, 
with surviving cool gas in the form of dense clumps with 
f b ~ 0.01 - 0.03. 

As is common elsewhere in astrophysics, the approach 
adopted in this paper is phenomenological modelling. A 
simplified picture of the phenomenon is combined with ap- 
proximate treatments of the expected physical effects to 
create an 'end-to-end' tractable theory that obeys conser- 
vation laws and makes testable predictions. Such theories 
are of course always an interim measure, to be discarded 
when the obstacles to calculation from first principles are 
overcome. Unfortunately, in this case, these obstacles arc 
formidable: 3-D time-dependent gas dynamics, radiative 
transfer, and heat conduction including saturation and pos- 
sibly magnetic suppression. 

Evidently, the fundamental approach is unlikely to yield 
results anytime soon. Accordingly, possible improvements 
of the crude modelling described herein should be inves- 
tigated. Also diagnostic codes should incorporate features 
of such models to extract more reliable parameters from 
observational data. 



Appendix A: Heat conduction into a spherical blob 

In zone 2, the blobs are surrounded by gas whose temper- 
ature is rising to coronal values. Conduction will therefore 
transfer heat into the blobs, and this constitutes a loss term 
in the energy equation for the ambient gas. 

Given that /& <C 1, it suffices to consider a single spher- 
ical blob with temperature Xf, and radius a located at r = 
in an infinite medium with T — > T a as r — > oo. If C is the 
cooling rate per unit volume and k is the conductivity, the 
equilibrium temperature profile for r > a is given by the 
equations 

^^4^[C(T)-C(T a )} (A.1) 

and 



dr ATtr 2 n 
with boundary conditions 

T(a) = T b and T(oo) = T a (A.3) 

In the blob's absence, the gas is isothermal and has cooling 
rate C(T a ), which is subtracted in Eq. (A.l). Accordingly, 
£(oo) is the additional cooling due to the blob's presence. 
Of this, AC = C(oo) — C{o~) represents emission from am- 
bient gas cooled below T a by the blob, and C(a) is the rate 
of heat conduction into the blob. 

In thermal equilibrium, C{a) is balanced by emission 
from within the blob. Now, since radiative cooling is cx p 2 
and pi, 3> p a , we expect that AC -C C(a). Therefore, to a 
first approximation, C(r > a) — C(a), a constant, and this 
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allows Eq. (A. 2) to be solved analytically when k oc T 5 / 2 
(Spitzer 1962). The resulting temperature profile is given 
by 

* 7/2 =< /2 + (l-< /2 )(l-^) (A.4) 

where t = T(r)/T a , and the corresponding rate of heat 
conduction into the blob is 

£d = ® Tra [ (kT)„ - (kT) 6 ] (A.5) 

Since kT oc T 7 / 2 , £ c ; is insensitive to Tb when T a ^> Tf,. For 
the solutions reported in Sects. 5 and 6, we take k = 1.0 x 
10 -6 T 5 / 2 , corresponding to Coulomb logarithm In A = 17. 

The above discussion treats conduction in the diffusion 
limit - i.e., where the mean free path of the electrons is 
<C macroscopic length scales. In the opposite limit, heat 
conduction into the blob is flux-limited and saturates at 

C sat = 47TCT 2 q sat (A. 6) 

where q sat is estimated by Cowie & McKee (1977) to be 

/2kT \ 1/2 
q sat = 0.4 n e kT e (A.7) 

V 7rm e / 

and is here evaluated at T a , (n e ) a . Interpolating between 
these limits (cf. Balbus & McKee 1982), we take the rate 
of heat conduction into the blob to be C con d, where 

^ = ^+^7at (A.8) 

At high temperatures, this gives L con( i oc T 3 / 2 in place of 
C cl oc T 7 / 2 . 

The above formula is for a non-magnetized plasma. But 
since stars in and near the weak-wind domain have detected 
magnetic fields (e.g., Oskinova et al. 2011), we include the 
possibility of magnetic suppression of heat conduction by 
writing 

Cm = 4> Ccond (A. 9) 

In this investigation, <fi is varied to explore its impact on the 
solutions. In future, it may be determined or constrained 
by fitting observational data. 

The suppression of thermal conductivity in astrophysi- 
cal plasmas has been strikingly confirmed by the discovery 
of cold fronts in X-ray maps of clusters of galaxies (e.g., 
Carilli & Taylor 2002). For the cluster Abell 2142, Ettori & 
Fabian (2000) estimate a reduction factor of between 250 
and 2500. They speculate that, as a result of merging, dif- 
ferent magnetic structures are in contact and so remain to 
high degree thermally isolated. The displacements of wind 
clumps from their nascent ambient surroundings might well 
lead similarly to substantial reduction factors. 
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